Lightly doped dimerized spin chain in the one-dimensional t-J-J' model 
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We study the one-dimensional t- J- J'-model in the limit of small hole doping x and small J/t, 
J' It. Special emphasis is put on the regime J'/ J « .5 where a spin gap is present at small doping 
and the undoped spin chain is strongly dimerized. Using a perturbative approach and Luttinger 
liquid arguments, we demonstrate for this non-integrable class of models that the charge degrees of 
freedom behave as non-interacting spinless solitons in the dilute hole limit. Our approach is also 
used to evaluate the energy and mass r enor malization of a single hole. Interestingly, the corrections 
of these quantities are in powers of y/j/t. At J' / J = .5 we construct a variational spin-polaron 
wave function for the hole and find good agreement with our perturbative results. 

PACS numbers: 74.72.-h, 75.10.JM 
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I. INTRODUCTION 

The discovery of high-T c superconductors has generated 
great interest in models of strongly correlated insulators 
in low dimensions, both in the presence and absence of 
doped carriers. While in two dimensions the problem 
remains challenging even for the most elementary mod- 
els that are believed to capture some of the physics of 
the cuprates, much can be learned from the study of 
one-dimensional counterparts of such models. Basic fea- 
tures that characterize the high-T c materials, such as the 
absence of a quasiparticle pole near the Fermi surface 
and possibly spin-charge separation, are essential to one- 
dimensional systems and are well understood herei. The 
success of the theory in one dimension is due to the avail- 
ability of powerful exact methods, such as numerical di- 
agonalization and Bethe ansatz, in combination with the 
knowledge of a low-energy effective theory which has so 
far been found to describe all gapless degrees of freedom 
in one-dimensional systems, the Luttinger liquid 2 . 

In the present work we study the effect of doping a 
small concentration of holes into a dimerized spin chain. 
This is motivated by the idea that in a dimerized system, 
there naturally exists an amplitude for pair formation. 
Upon introduction of carriers phase coherence may be es- 
tablished, resulting in superconductivity. This idea has 
been previously proposed in the literature^. We distin- 
guish between two rough physical pictures of hole doping 
into a dimerized chain (Fig. In the first scenario (Fig. 
I^i)) the dimer order remains long ranged even after hole 
doping. The dimers tend to reside on every second link 
of the lattice as in the symmetry broken undoped state, 
and holes tend to pair on the empty sites in between. 
We note that this case bears some resemblance to doped 
ladder models^. Indeed strong superconducting correla- 
tions have been proposed for dimer-models with explicit 
symmetry breaking^. 

A second possibility is that the holes enter as domain 
walls between different dimer phases (Fig[IJ)). Then the 
long-range order of dimers is destroyed, but the singlet 
gap remains. If the holes are mobile and become phase 
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FIG. 1: a) doping into frozen dimer state. Ovals represent 
singlet pairs, b) mobile dimers with domain walls. 

coherent, this leads to superconductivity in analogy with 
the original RVB scenario^. However, the pictures pre- 
sented in Fig. are based on the limit of large J/t. In 
physical systems, on the other hand, J/t is usually a frac- 
tion of unity. In this regime, the hole kinetic energy is 
dominant and it is less clear intuitively if the dimer pic- 
ture discussed above applies. In the present paper, we 
study this question in the context of the t- J- J'-model in 
one dimension: 

H = -tPy^(ct .c, + i |g + h.c,)P + J y^QS, • S i+ i - -mni+i) 

i i 

+ J f y](Sj • S i+2 - ^nin i+2 ) + y y^n i+i 

i i 

(i) 

Here, P is a projection operator that enforces the con- 
straint of no doubly occupied sites. The Si are spin-1/2 
operators, and the nearest- and next-nearest-neighbor 
couplings J and J' are assumed positive throughout this 
paper. In addition, we have included a nearest neighbor 
interaction V for later convenience. 

At finite doping, not much is known analytically about 
this model due to its non-integrable nature. The phase 
diagram has been established numerically^^ (Fig. EJ). 
However, in the regime around a = J' / J ~ .5 where the 
undoped spin chain is strongly dimerized and which is of 
particular interest to us, we feel that the numerics are 
somewhat inconclusive for small doping x and J/t < 1 
as we will discuss below. 
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FIG. 2: Sketch of the zero temperature phase diagram of the t- J- J' -model as determined numerically in Ref0 for a — (a), 
a = .24 (b) and a = .5 (c). Contours are labeled by values of K p . The shaded region marks the domain of dominant singlet 
superconducting correlations. The dotted lines in (a) were proposed in Ref. 0. 



The purpose of this paper is twofold: First, in the case 
of strong frustration a = J' / J sa .5 we wish to deter- 
mine the fate of the various regions present in the phase 
diagram of the t-J-J'-model at small doping and small 
J/t (FigEt) ■ In particular, we will answer the question 
whether a regime of dominant singlet superconducting 
correlations in the vicinity of an instability towards phase 
separation persists to values J/t <C 1 and in the dilute 
hole limit. This would happen if all the contours in Fig. 

extrapolate to the origin, which appears to be a pos- 
sible interpretation of the numerics. To this end, we in- 
troduce a perturbative approach valid in this limit. We 
find that the scenario mentioned above does not occur, 
but instead the Luttinger liquid is stable for small J/t 
and small doping. Luttinger liquid arguments will then 
imply that strong superconducting correlations only ex- 
ist above a finite critical value of J/t. 

Second, below the critical value (J/t) c the liquid phase 
is stable in the limit x — ► 0, and we will use our method 
to demonstrate certain properties that one expects to 
hold for one-dimensional lattice models based on general 
grounds. In particular, spin and charge are expected to 
correspond to separate degrees of freedom and any micro- 
scopic coupling between them should be irrelevant. As a 
consequence, dilute holes that are doped into a correlated 
spin chain should act as a gas of non-interacting spin- 
less solitonsii, where the coupling to the non-trivial spin 
background only gives rise to a renormalization of the 
effective hole mass. This phenomenon has been observed 
in integrable models such as the Hubbard models. Our 
perturbative approach allows us to give a demonstration 
of the same behavior in the non- integrable t-J-J'-model, 
and to calculate the effective energy and mass renormal- 
ization of the hole for small J / 1. 

The remainder of this paper is organized as follows: 
In section ITT1 we introduce the model, briefly discuss the 
numerical phase diagram and cast the model into a lan- 
guage where holes are interpreted as domain walls. In 
section ITTT1 we treat the spin-charge couplings as a per- 
turbation and derive expressions for the ground-state en- 



ergy, compressibility and Kohn stiffness of the model for 
small J/t, J'/t and small doping x. This will allow us to 
qualitatively continue the numerical phase diagram into 
the region of small doping. In section IIVI we will ex- 
plicitly evaluate these expressions as asymptotic series in 
powers of x and */ J/t, demonstrating the convergence 
of our approach. In section E] we discuss the single po- 
laron energy and mass renormalization. We introduce 
a spin polaron type variational wave function where the 
hole is surrounded by a cloud of tightly bound triplet ex- 
citations and find good agreement with our perturbative 
results. We conclude in section IVTl Appendix ^] illus- 
trates the behavior of our expansion at general order in 
various limits. An important technical issue is discussed 
in appendix 151 



II. FORMULATION OF THE PROBLEM 

We wish to study the Hamiltonian Q in the limit of 
vanishing exchange couplings and doping. For V = 
the phase zero temperature diagram of Q has been ob- 
tained numerically^iiS for various values of the parame- 
ter a = J' / J. The results are sketched in Fig. These 
phase diagrams show a Tomonaga-Luttinger liquid (TL) 
region below the dashed line labeled spin gap, where both 
spin and charge degrees of freedom are gapless. Above 
the dashed line there is a spin gapped liquid phase which 
is subdivided into a regime of dominant singlet super- 
conducting (SS) correlations (shaded grey) and, where 
present, a regime of dominant charge-density- wave cor- 
relations (CDW). Also shown are contours of constant 
values of the Luttinger parameter K p . Above the K p 
line labeled oo, there is a region of phase separation (PS), 
where the ground state has a phase boundary between a 
hole-rich and an electron-rich phase. The parameter K p 
is directly related to the large distance behavior of the 
various correlation functions of a Luttinger liquid: For 
K p > 1 pairing correlations dominate over density-wave 
correlations, otherwise density- wave correlations domi- 
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nate. Furthermore, in the presence of a spin gap, triplet 
pairing correlations and spin-density-wave correlations 
are exponentially suppressed. In this region, the main 
competition is therefore between CDW and SS correla- 
tions, which in the presence of a spin gap decay as r~ Ap 
and r~ 1 / Kp respectiveljii^. 

At zero doping (x = 0) it is well known that the spin 
chain undergoes a phase transition at a critical value of 
m .24 (Ref. [H above which the ground state is dimer- 
ized with a gap in the spin excitation spectrum^. For 
a > a c this spin gap remains present over a finite range 
of doping for any value of J/t, as was shown in ReflloL 
Fig. [3 shows that at a = .5 the spin gapped region has 
considerable overlap with the pairing region K p > 1 even 
at small values of doping x. Furthermore, it appears from 
the numerics as if all K p contours, including the phase 
separation boundary K„ = oo, flow to small values of J/t 
at small x (cf. Refs. I4 ¥l0r) . One expects that these con- 
tours will focus on a critical point [x = 0, (J/t) c ], as was 
proposed in Ref. for the t- J-model (J' = 0, Fig. |2K))- 
A possibility that seems consistent with the numerics at 
a = .5 is that (J/t) c = 0, i.e. all contours flow into the 
origin of the phase diagram. In this case, a sufficiently 
small amount of doping would always lead to phase sepa- 
ration, and upon further doping one would enter a region 
of dominant SS correlations. Alternatively, (J/t) c could 
be finite but possibly smaller than its value at a = 0, 
which is between 3 and 4 (Fig. This would im- 

ply that the above phenomenology of phase separation 
and superconductivity at small doping occurs only for 
J/t > {J/t) c , while for J/t < (J/t) c the liquid phase is 
stable at any doping x. In the latter case, one would ex- 
pect the Luttinger parameter K p to approach the value 
1/2 in the dilute hole limit x — > 0, which is the value 
corresponding to non-interacting spinless degrees of free- 
dom. This behavior is clearly exemplified by the nu- 
merical phase diagram of the t- J-model (Fig. [2^). The 
primary goal of this paper is to determine which of these 
two scenarios applies to the t-J-J' model at a > a c . 

We begin our analysis by casting the Hamiltonian 
into a language where the holes play the role of domain 
walls between broken segments of an infinite spin chainSS 
(Fig. |2J. We consider a lattice of L sites with a number 
of N e electrons and Nh = L — N e holes. Denoting the 



i'th spin on the lattice by Si we may regard the spins as 
residing on a "squeezed" lattice where the hole sites have 
been dropped from the system and the label i of the spin 
Si is a site label in this squeezed space, as in Fig. |3J). We 
also introduce interstitial sites for the squeezed spin lat- 
tice whose labels j = i + 1/2 differ from those of the spin 
sites by 1/2. Each interstitial site may accommodate a 
number rij = 0, 1... of holes. A faithful representation of 
the Hilbert space of is given by states labeled by 

• ■ .,(Ti,n i+ i,a i+ i, ...^ = \...ai...)\...nj...) 

(2) 

i =l...N e 



where <Ji = ±i denotes the z-component of the spin S%. 
This language turns out to be particularly convenient 
when one introduces a large nearest neighbor hole repul- 
sion, i. e. if we let V = oo in 0J, such that it is forbidden 
for two holes to occupy neighboring sites. In the present 
language this means that the occupancy of the hole clus- 
ter labeled by j is now restricted to be rij = 0, 1. This 
modification of the model will be irrelevant in the dilute 
hole limit which we are interested in2£. We may now 
choose to formulate the hole dynamics either in terms of 
hard-core boson operators or spinless fermion operators. 
For convenience, we introduce fermion ladder operators 
Cj , Cj , where the action of cj can be described by cutting 
the spin chain open at the interstitial site j, introducing 
a hole at this position and multiplying the state by an 
appropriate fermion phase. The hole kinetic energy is 
then simply given by 

H c = -^2(tc]cj +1 + h.c?) (3) 

3 

© can be thought of as the J = J' = limit of the 
Hamiltonian. In the other limit of interest, namely the 
limit of zero doping x — Nh/L, the Hamiltonian becomes 
that of a pure spin chain: 



i i+1 i+2 

x x O x x x O 

j=i+l/2 j+1 j+2 

FIG. 3: Domain wall representation of the Hilbert space. The 
upper chain represents the spin sector with sites labeled by 
i. The holes live on a lattice of interstitial sites labeled by j. 
Crosses represent "empty" interstitial sites, circles represent 
holes between two spins. 



H S = J^2 {Xi,i+i + aX hl+2 ) 

1 (4) 

where X i V = Si - Sy - - 

4 

where we work at constant a = J' / J from now on, and 
assume that a > a c , such that the small doping regime 
is spin gapped. The combined Hamiltonian 

H =H S + H C , (5) 

where the spin and charge part are still completely decou- 
pled, will serve as a starting point for the perturbation 
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theory we propose. In order to correctly reproduce ma- 
trix elements of the Hamiltonian ft|. couplings between 
the spin sector and the charge sector must be introduced: 



H = H + H ac + H' sc (6) 
H sc = -J^JVYj (7) 

3 

H' sc = Ja^Xi-^i+w^iUi+i (8) 

i 

where 

7 i+ i = (1 - a)^M+i + + A M+2 ) (9) 

Here, i? sc is a correction which couples spin and charge 
by adjusting nearest neighbor bonds and removing next- 
nearest neighbor bonds in the squeezed spin space in the 
vicinity of a hole. Certain corrections of the latter sort 
are redundant whenever two holes are next-nearest neigh- 
bors in real space - or nearest neighbors in the present 
formalism - and this is corrected by H' sc . Formally, H sc 
and H' sc are suppressed by powers of both J/t and x, and 
hence can be regarded as small compared to Ho- Our 
strategy is thus to treat the spin-charge coupling terms 
H sc + H' sc as a perturbation. We must caution, however, 
that the hole kinetic energy H c is very small, or order 
tx 2 , and the small energy denominators that appear in 
perturbation theory must be treated with care. We apply 
this method to the spin gapped regime a > a c and find 
it to be a valid procedure in second order perturbation 
theory, in the sense that corrections are indeed small in 
the limit we consider. The general systematics of this 
at higher order perturbation theory are elucidated in ap- 
pendix ^ 

In the following section we will show that our approach 
gives rise to a perturbative expansion of the ground-state 
energy which may be used to analyze the phase diagram 
of in the vicinity of the origin, where the fate of the 
various phases is uncertain from numerics for a > a c 
(Fig'Ht). We note that the procedure proposed here bears 
some resemblance to that used by Xiang et al. to study 
the t- J-model in a first order perturbative approach^. 



III. PERTURBATIVE ANALYSIS OF THE 
MODEL 

In gapless one-dimensional systems it is generally possi- 
ble to derive basic features of the phase diagram from 
spectral properties by means of Luttinger liquid theory 2 . 
The low-energy properties of a Luttinger liquid are com- 
pletely defined in terms of three parameters which have 
the dimension of a velocity: A sound velocity v s , a "com- 
pressibility" parameter wjv related to elementary charge 
excitations, and a "stiffness" parameter vj related to ele- 
mentary current excitations. These are not independent, 
but are related by the following universal relations iden- 



tified by Haldane^: 



vn = v s /Kp, vj = v s K p 



(10) 



The parameters vn and vj can be calculated from the 
dependence of the ground-state energy Eq on the carrier 
density x and on a phase twist <j> respectively^, where <fi is 
associated with a flux Lip penetrating the system when it 
is imposed on a ring with periodic boundary conditions: 



vj 



2 d 2 E 
irL dx 2 
7T d 2 E n 



2L d<j) 2 



(11) 
(12) 



6=0 



By (lllfl . (|12fl . vn is proportional to the inverse com- 
pressibility of the system, while vj is proportional to the 
Kohn stiffnessii, which is related to the Drude weight of 
the conductivity. Eqs. (|1U|) - (|12H allow the determination 



of K p via 



(13) 



from the ground-state energy of the system alone. The 
strategy in now to evaluate both the numerator and the 
denominator in l|13(l perturbatively. 

We now proceed by imposing periodic boundary con- 
ditions on the charge sector and the spin sector of the 
system (0-© separately. This is apparently not the 
same as imposing periodic boundary conditions in real 
space, since momenta are now quantized in units of 
2ir/N e rather than 2n/L. Note that there is a unique 
and well defined map between the state space introduced 
in (01 and the Hilbert space of the t-J-J' model only for 
a finite system with open boundary conditions. Indeed, 
imposing open boundary in real space is equivalent to 
imposing them in the spin sector and the charge sector 
separately. However, going from open to periodic bound- 
ary conditions is not expected to matter for large system 
sizes. The unperturbed Hamiltonian Hq then has two 
separately conserved momenta, and we denote its ground 
state by 



\a a ,ip ) = Wo) |"0o> 



(14) 



where | cr ) is the ground state of the spin Hamiltonian 
H s on a ring of N e spin sites. Although for a > a c the 
ground state of H s has a broken translational symme- 
try and is doubly degenerate, we will assume that |ao) 
is a symmetric superposition of the two symmetry bro- 
ken ground states and thus has zero lattice momentum. 
Likewise |"0o) is a non-interacting Fermi sea of 2V/, spin- 
less Fermions hopping on N e sites with periodic bound- 
ary conditions. The unperturbed ground-state energy we 
write as 
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2t 

E =E ao +E i)a =E <T0 -N e —sm(k f ) (15) 

7T 



N h TTX 



where kf = ir-rr — 



N P l-x 



(16) 



and £"0-0 , are the ground-state energy of H s and the 
spinless fermion kinetic energy, respectively. We will fo- 
cus our analysis on the limit x — > where J/t is small 
but fixed. In this limit we argue that the ground-state 
energy of © has an asymptotic expansion of the form 



E = E aa + L(Ax + Bx 2 + Cx 3 + ...) 



(17) 



The coefficients A,B... will depend on J/t. At the lead- 
ing order, they can in principle be inferred from the spin- 
less fermion kinetic energy in JTSJ. Formally, however it 
will be more convenient to work with an expansion of the 
form 



E = E, 




The coefficients in (|17|) and <|18ll will in general not be 
the same, due to the non-linear dependence of kf on a; in 
(|16fl . However, since N e = L(l — x), the A term in i|18|) 
is linear in x, and hence 



A = A, B = B 



(19) 



We will now proceed by evaluating the above series order 
by order in perturbation theory, treating the spin-charge 
coupling terms H sc + H' sc perturbatively as we have out- 
lined above. We write 



E = E +E 1 +E 2 + ... + E k + . 
A = A Q + Ai + A 2 + . . . + A k + 



was used. Eq. (|22|l shows that if B acquires a finite 
negative value due to the spin-charge couplings, K p will 
diverge as x — > even at small J/t. This would imply 
strong superconducting fluctuations and phase separa- 
tion at the divergence. On the other hand, if B is zero 
or positive, the liquid phase will be stable for small x 
and J/t, and dominant superconducting correlations will 
be absent in the vicinity of the origin of the phase di- 
agram. In this case, one would expect that B = to 
all orders in perturbation theory, since for B > one 
would have K p — > as x — > which seems inconsistent 
with the numerical phase diagram. Also, K p = is a 
somewhat unlikely pathological limit of Luttinger liquid 
theory, where the coefficient of the conjugate momentum 
of the charge field vanishes and the charges freeze into a 
classical state. This again seems unlikely in the absence 
of long-rage interactions, which one may assume espe- 
cially in the presence of a spin gap. As we have argued 
initially, one would rather expect the charges to behave as 
non-interacting spinless solitons in the dilute limit, and 
Luttinger liquid physics then implies that K p assumes 
the value h in this limit, provided that no instability in- 
tervenes. This then requires that the coefficient A^ and 
C in \T2l are not independent, but have a constant ratio 
independent of J. 

Hence, we distinguish only two cases which we feel are 
the only ones consistent with the numerical phase dia- 
gram for a > a c , Luttinger liquid theory and general 
expectations from the study of integrable model o 2 i 12 : 1. 
B < 0, leading to a phase separation instability in the 
dilute hole limit for any value of J/t, and 2. B = 0, 
6C = tt 2 A^ for any J, corresponding to a stable liquid 
as x — > and small J/t, where the charges act as a di- 
lute gas of non-interacting spinless solitons. We will now 
show that the second case applies. 

At first order perturbation theory the energy correc- 
tions factorize into mean-field like products, since spin 
and charge are not correlated in the ground state wave 
function 114(1. We have 



(20) 



and similarly for all other coefficients, where the label k 
denotes a term arising at /c'th order perturbation theory. 
We have: 



hence: 



Ei = (cr , 4>o \H SC + H' sc \ cr , V'o) 
= —N e J (rij} (7j) + N e Ja (X f 
= -L(jj) Jx + LJO(x 4 ) 

Ai = - (7i> J 



-M+i) ( n i n j+i)o 
(23) 
(24) 



A Q = -2t B = C = —t (21) 
From (CTJ- (USJ , 113 it follows that at small J/t 

7T / A^X _ 1 / tx 

p ~ 2\2B + 6Cx ~ 2y B/^+tx [ > 

where A^ denotes the second derivative with respect to 
the phase <fi introduced above, and Aq(c/)) = —2tcos((f>) 



where <>o denotes the expectation value with respect to 
|(To > or IV'o > when no ambiguity is possible. Note that 
the contribution of H' sc is of order a; 4 . The smallness of 
this term as x — *■ reflects the fact that the holes obey the 
Pauli principle which suppresses the probability of two 
holes being near each other. We see that already at this 
order, H sc does not renormalize any of the coefficients 
in l|17|l that we are interested in. We will thus drop it 
from the subsequent discussion. To determine leading 
corrections to B and C, we will need to go to second 
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order: 



terms into a double integral over a hole momentum k\ 
and a particle momentum fc 2 : 



En 



(vo^ol H sc W,tf}) W,tp\ H sc \a ,tpo) ^ 



Eif, — E^ + E„ — E a 



where the sum goes over a complete set of unperturbed 
eigenstates and the prime excludes the ground state lfH|) 
from the sum. We now rewrite H sc as 



H, 



— V 

N r ^ 



n q j- q 



(26) 



E' 2 = -N £ 



dki 
~2n 



2tt 



f(h,k 2 ) (31) 



where 

i ct ; 

and e(fc) = — 2tcos(fc) is the free fermion dispersion. For 
later convenience, we also introduce the function 



where Fourier transforms 
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E 



(27) 



2 C^2' 



J 2 
2A^ 



(o-ol7- 



l_7g Icq) 



-Evt — Er, 



+ {<i 



(33) 



have been introduced. Using the fact that the intermedi- 
ate states in (|25|1 can be chosen to be momentum eigen- 
states, we have 



where the symmetry of e(k) was used. 

The leading correction to the energy at second order 
perturbation theory is a contribution to the A coefficient 
in jnj|: 



Ei 



(ipo\n- q \ 

^2 2.^1 2.^1 

e i \<y,i>) 



n q \ipo) (cr |7? W) (c\l-q Wo) 



Ej, — E, 



Er, — Er, 



(28) 

It is necessary to distinguish between terms with zero 
momentum exchange between spin and charge and those 
with q 0. We write: 



En 



En + En 



(29) 



where E® contains all q — terms and E' 2 contains all 
the rest. At q — 0, n q — J2 n j ~ commutes with 
the Hamiltonian, hence there can be no virtual charge 
excitation and the charge matrix element is diagonal, 
\tp >= |V>q >: 



El 



J 2 ^ 



(g |7g^o W) (g1_7g=o Wo) 

Err — E ar . 



(30) 



Note that virtual states without spin excitations do not 
enter Q28[l. since \a >= Wo > would imply q = and 
again the charge part vanishes unless also \ip >= \i/}q >, 
which is excluded from the sum. Thus for a > a c the 
energy denominator in l|28|l is bounded from below by the 
spin gap A which will dominate over charge excitation 
energies of order tx 2 very close to the Fermi surface. This 
assures that the perturbative expansion is well behaved 
in the limit x — > (see appendix 0). 

For q we note that n q excites only single particle- 
hole excitations. We can thus convert the sum over these 



k d , 

A2= W e dk~ f E2 



2jv 



dk 

2^ 



(34) 



f(0,k) 



To leading order in 3 ft the integral over momenta may be 
carried out to give a quantity defined in terms of the pure 
spin chain H s . This will demonstrate that the present 
expansion is well behaved, but we defer the evaluation 
to the following section in order to continue with the 
analysis of the crucial B coefficient. Its correction at this 
order reads: 



Bo 




-F(0) 



2N e dk 2 2 



(35) 



k f =0 



Again, the contribution from E' 2 is evaluated by straight- 
forward differentiation of l|31(l . Only boundary terms sur- 
vive, as all derivatives of the integrand vanish by sym- 
metry when the limit kf — > is taken, using the 27T- 
periodicity of / in the second argument. We find: 



B 2 = -F(Q)+ \(F{2k f )+F(2ri)) 



(36) 



At this point we have introduced an infinitesimal rj since 
terms with zero momentum transfer are really excluded 
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in the sum defining E' 2 . However, we argue that the 
function F(q) will be continuous at q = and hence B 2 
vanishes. We note that this is the effect of a non-trivial 
cancellation between q = processes and processes with 
g ^ 0. Physically, the continuity of F(q) can be seen by 
interpreting F(q) as the second order energy response of 
a pure spin chain due to a periodic perturbation, as we 
explain in detail in appendix IB! Hence, by (|36|) 



following replacement of the free hole dispersion in the 
function /(fci,^) in 1|32[1 : 



e(k) — ► e^(fc) = — 2<cos(fc + . 
f(ki,k 2 ) — ► f l p(k l ,k 2 ) 



Hence from (|34|) 



(41) 



B? = 



(37) 



We have convinced ourselves that a similar cancellation 
in the B coefficient takes place at third order perturba- 
tion theory^. We therefore propose that 



B = 



(38) 



to all orders in perturbation theory, and thus for small 
J/t the liquid remains stable in the limit x — > even 
in the case a ss .5 (Fig. 0:). The physical implication 
of l|38|) is that indeed the holes act as spinless fermions 
whose interaction is short ranged, and is irrelevant in the 
dilute limit. The Pauli principle severely suppresses the 
wave function when two holes approach each other. The 
range of this suppression is larger in one dimension than 
for dimensions greater than one, since in higher dimen- 
sions a curvature of the wave function is less costly at 
small distances. Therefore, in one dimension this effect 
is strong enough in order to prevent a short range inter- 
action from generating a term of order x 2 in the energy. 

As we have argued above, the non-interacting nature of 
the charge degrees of freedom in combination with Lut- 
tinger liquid arguments also imposes constraints on the 
linear and cubic terms in x when a magnetic flux is im- 
posed. We now move on to verify these relations pertur- 
batively. Note that B = leads to C = C in (JT7I) and 
(jT%|l . such that 



^2,0,0 — Tr7T^2 



I s W MQ - k) 



However, using the fact that 



',=0 

(42) 



(43) 



holds, it follows by shifting the integration variable and 
comparison with (gOJ that 



6C ? = ir 2 An 



(44) 



is satisfied. Again, we have confirmed an identical rela- 
tion at third order perturbation theory^, and this sug- 
gests that indeed 



6C = tt 



(45) 



to all orders. Hence, although the parameters u/v and 
Vj in l|13() each receive nontrivial corrections, their 
ratio is fixed to leading order in x such that K p always 
approaches \ in the limit x — > 0. Luttinger liquid theory 
then implies that the dilute holes share all the universal 
properties of a gas of non-interacting spinless particles. 
Presumably, this picture will hold in the entire regime 
J/t < (J/t)c, where for J/t > {J/t) c small doping will 
give rise to phase separation. 



C 2 = 



6N e dkj 2 



(39) 



fe f =0 



In this case boundary terms such as those displayed in 
H36JI do not contribute, since they vanish by symmetry 
as kf — » when another derivative is acting on them. 
Instead, we have now again a "bulk" contribution analo- 
gous to that in (|34|l : 



C 2 



1 7T 

6 2 



271 



(40) 



k f =0 



As explained above, this is to be compared to the coef- 
ficient ^2.00. The phase twist <p w iU modify all hopping 
matrix elements via t — > t e 1 ^ in © and leads to the 



IV. EXPLICIT EVALUATION OF 
COEFFICIENTS 

In the preceding section we have shown that our pertur- 
bative approach is consistent in all details with a pic- 
ture where the charge degrees of freedom behave as non- 
interacting spinless solitons in the dilute limit, and are ef- 
fectively decoupled from the spin dynamics. The second 
order expressions we derived involve complicated sums 
over both spin and charge degrees of freedom. We will 
now show that the expressions for A 2 and C 2 can be eval- 
uated more explicitly, to leading order in J/t, in terms of 
quantities that are derived from a pure spin chain prob- 
lem. In this way we obtain explicit asymptotic expan- 
sions for the ground-state energy and the compressibility 
parameter ujy, showing that second order corrections are 
suppressed by non-trivial powers of J/t compared to the 
leading orders. Also, these quantities are related to the 
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single hole energy and mass renormalization, which will 
be clarified in the following section. 

We stress once more that the results we present here 
are valid in the limit x 2 -C J/t. In this limit the low 
lying charge excitations are dominated by the curvature 
near the band bottom of the bare dispersion e(k) and 
their contribution to the energy denominator in l|31fl is 
dominated by that of the gapped spin excitations. In 
the opposite limit J/t <g; x 2 the perturbation theory pre- 
sented here is still valid, yet a crossover will take place 
and the asymptotic expansion (|17fl will not hold (see ap- 
pendix . 

With this in mind, the first and second order energy 
corrections are dominated by the following terms: 



Ex/L 
E 2 /L 



A 2 x 



Jx 



(46) 



We will now show that the second order term is indeed 
suppressed by powers of J/t compared to the first or- 
der term, which is of order Jx. To achieve a systematic 
expansion of A 2 in J/t we rewrite <|34|) in the form 



~ J dEj2K n (E) 



(51) 



where the coefficient Ti is a quantity defined only in 
terms of eigenstates of the doped spin chain. For later 
convenience we define the more general function: 



r P = J p ' 2 



dE 



A(q = 0,E) 
{E-E ao f 



(cTo|7«=ok)<a| 7g =oko) (52) 



E„-E„ n \P 



Hence, it is apparent from l|46(l and (|51J) that 



E 2 
E~i 



(53) 



A, = -l "£/,/£ 



A(k,E) 



e(k) - e(0) + E - E a 



where we have introduced a spectral function 



A(k,E) = w J2 \(^\lk\a)\ 2 S(E a -E) 



Y^K n {E)e 



(47) 



(48) 



and its energy dependent Fourier coefficients K n (E). In 
terms of the latter we may write 



indicating the convergence of our perturbative approach 
for small J/t. In appendix ^ we will further comment 
on convergence and expansion parameters of this series in 
various limits. Note that the non-analytic nature of the 
expansion originates from the gaplessness of the charge 
degrees of freedom and the existence of a regime where 
the spin gap dominates the energy denominator in 1471 . 
Similarly, in second order perturbation theory the com- 
pressibility parameter reads to leading order in x: 



v N ~ ( Ant + — C 2 ) x 

7T 



(54) 



By means of (f40|) the evaluation of C 2 goes analogous to 
that of A 2 and we get: 



+ - 9 ^ g ° — cos(&) 
(49) 



The fc-integral is readily performed to give 



A 2 - 

where A = 



1 f „ (l + A-V(l + A)2-l 
- / dE^KniE)^ 

E -E, 



V(l + A) 2 -1 



2t 



(50) 



The matrix elements defining K n {E) will decay rapidly 
when E is a few times J and hence we may expand l|5U|) 
in powers of A. Keeping only the leading term this yields 



v N ~wt u-rw- \x 



(55) 



Numerically, we found Ti = .2502(2) and T| = .474(3) 
at a = .5 (Fig. @J. Hence, although the compressibility 
k ~ v^ 1 increases with J, no unstable value of J can be 
inferred that lies within the validity of our perturbation 
theory. 



V. SINGLE SPIN-POLARON PICTURE 

We will now develop a variational picture of the polaronic 
effects of a single hole on its spin environment at small 
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p=0 
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0.1 




average 









fit 





0.05 ' 

0.05 



0.1 1/N 




0.1 1/N 



FIG. 4: Numerical determination of Ti and F3 by exact di- 

2 2 J 

agonalization of H s for a = i. System sizes of up to N = 18 
have been diagonalized. Results are plotted for the two de- 
generate ground states with momenta p — (crosses) and 
p — 7r (squares). The extrapolated values have been deter- 
mined by fitting their averages (diamonds) to the function 
f(N) =a + b exp(-cN) 

J/t in the special case a — .5. The perturbation the- 
ory presented in the preceding section for a finite carrier 
concentration may be applied to the problem of doping a 
single hole into the infinite spin chain as well, such that 
we will be able to compare variational and perturbative 
results. In second order perturbation theory, the energy 
of a single hole at momentum k reads: 



E p (k) = -2tco S (k)-( 7j ) J- I f(k,k 2 ) + 



(56) 



where we have not included the contribution E ao from 
the spin background. At k = we immediately see by 
comparison with (|2*T)l . (|2"4"|) and that 



E p = E p (k = 0) = A + A 1 +A 2 + ... = A (57) 



holds for the single polaron energy in second order per- 
turbation theory. Likewise, for the renormalized mass of 
the spin polaron we have, comparing to eqs. 1)21(1 and 
60)1 



m 



dk* p 



k=0 



— (C 



c 2 



)^^C (58) 



at this order. We may therefore rewrite the ground-state 
energy of the system at finite doping (|17fl as 



E = E n 



L ( E p x 



6m 



r k f fib 

E ao +N e J ^ 7^E P (k) +0{k 



(59) 



Hence up to third order in x the ground-state energy 
of the system is apparently given by the energy of non- 
interacting spinless particles with a dispersion E p (k), 
where interaction effects enter only beyond this order. 
This further confirms the picture established in the pre- 
ceding sections. 

We now focus on the Majumdar-Gosh point a = .5, 
where the ground state of the spin Hamiltonian H s is 
known exactlyiS. It consists of a direct product of un- 
corrected singlet pairs: 



i mg >=ii^(in>-ut) 



2i,2i+l 



(60) 



Note that we use \MG > to denote one of the two dou- 
bly degenerate symmetry broken ground states, whereas 
| Co > has been used to denote their symmetric superposi- 
tion. At a = .5, our results for the single polaron energy 
and mass eqs. (|57f) and l|58|l take the concrete form 



(61) 



We may write the ground state of the unperturbed 
Hamiltonian Hq as a superposition of states depicted in 
Fig. Et)+b): 




1 



(62) 



where \j) denotes a state with a hole at the interstitial 
site j and "c" and "s" refer to the spin sector and the 
charge sector of the state. When the interaction H sc is 
taken into account, a hole in the state \2j) will excite 
a spin configuration where the two dimers adjacent to 
the hole are in triplet states, and the two triplets form 
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a) 



b) 



c) 



d) 



2j 





2j 




FIG. 5: a)+b) Single hole basis states forming the ground 
state of the non- interacting Hamiltonian Ho (eq.J^J)- Lines 
denote the singlet pairs in (16011 . c) A pair of triplets excited by 
the presence of the hole. The oval denotes a singlet formed 
by two triplet states on the links adjacent to the hole, as 
displayed in 1631 . d) States used to form the variational wave 
function 164H . 

a singlet (Fig. HJ;)). More precisely, the oval in Fig. [SJ;) 
denotes the following spin state: 



(63) 



/,/ 



v ^ITT>« ®|U) W + IU>« ®ITt) w 
-^(|U> + I4T>)..® (|U> + I4T) 



We may now denote such a triplet pair excitation which 
is centered around the hole site 2j by \2j) . Similarly, 
the hole state \2j + 1) will excite the spin states \2j) 
and \2j + 2) s . Clearly, a variational wave function will 
need admixtures of states such as shown in Fig. 0:). 
However, in the vicinity of the hole the kinetic energy H c 
is the dominant part of the Hamiltonian, and it will allow 
the hole to move away from the excited triplet states as 
shown in Fig. 01). To optimize the kinetic energy, it is 
hence necessary to include the more general states of Fig. 
01) into the wave function. We therefore write down the 
following trial wave function for a state with one hole at 
momentum k: 



fa) 



\MG) S + e-l'-^-K 3 "-"') 4 \2f) 

2? 



(64) 



Hence \tk) lives in the subspace of all states that can 
be reached by acting once with the perturbation H sc on 
the zeroth order wave function l|62|) and then acting an 



arbitrary number of times with the hopping operator H c . 
We also note that indeed certain exact excited states of 
a Majumdar-Gosh spin chain are given in terms of the 
tightly bound triplet excitations shown in eq. (|63l) and 
Fig- El (Ref. \2(j) . The following matrix elements are 
needed to evaluate the energy of the state : 



,{MG\-J l3 \2]') 
S (2j'\-J 7j \2f) 



- ~£-J i S 3,2j' + S 3-l,2j> + tfj+W ) 

(65) 

= — J- (-1)3 — J 

~ 16 ( > 16 

(66) 



,(2j'|fl.|2i'),= 8 {MG\H S \MG) S + J 
= -^JN e + J 



(67) 



In addition, both and H s do not have off-diagonal 
matrix elements among the states \2j) s . This leads to 
the following expectation values: 



(t k \H c \t k ) = -2icos(fc) -4i|a| 2 £ 3 cos(fc + <5) 
(tk\H s \tk) 

(tk\H sc \tk) 



(68a) 



-JN e (t k \t k ) + J|a| 2 (ei+6) (68b) 



■■ ^J{t k \tk) 

— J (1 + 2e-P costf)) (a + a*) + 0( J\a\ 2 ) 
16 

(68c) 



(tk\t k ) = 1 + |a| 2 (a + 6) 



(68d) 



where the constants £i and £2 are proportional to the 
weight of spin excited states with the hole on even posi- 
tions and odd positions, respectively, and £3 arises from 
hopping between even and odd sites in the presence of a 
spin excitation: 



61 ^-Ve* 1 = — + -13 + 



a 5^ 4/3 6^ 



6 = \j2 e ~ m ' 1 (e-f , W- 1 \+e-f > W+ 1 \) 



4/3 24 A 
(69) 



Terms of order J\a\ 2 were only kept in (|68() when they 
are multiplied by £j ~ 1//3- It is apparent from l|68a|) 
that 5 = — k has to be chosen, and from l|68c|) that a is 
real and negative. Keeping only leading terms, this leads 
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to the variational energy function 



E v p ar (k; a, (3) = (t k \H c + H s + H sc \t k ) / (t k \t k ) 

9 1 a 2 

= -2t cos(fc) H J + -ta 2 [3 - t—(l - cos(fc)) 

w 16 2 ' P K v " 



+ J ^ + ^ Ja(1 + 2C ° S(fc)) 



(70) 



where again the bulk contribution of the spin chain was 
not included. We first minimize this function for k = 
and find for the variational parameters at the stationary 
point 



A> = \ T 



a 



J 

3\/3 17 



(71a) 
(71b) 




By ()71a[l . the size of the spin polaron cloud is propor- 
tional to (J/£)~ 5 in agreement with l|56l) where the dom- 
inant contributions to the integral come from the region 
where k\ is of the order of the spin gap. The variational 
energy of the spin polaron at k — is thus 



= £««•(*. = (); ao, A,) =-« 



This is indeed of the same form as (|61|1 where the first 
two terms are reproduced exactly, as they are mean-field 
like in character. Moreover, the coefficient of the last 
term is about .105 and hence matches the one obtained 
by perturbative and numerical methods in l|61|l within 

roughly 15%. 

The appearance of a mass term proportional to \J J/t 
as in (|61|l may also be understood from this variational 
approach. It is seen in 170JI that a term of order a 2 / '(3 ~ 
\f~J is no longer precisely canceled at finite k. The reason 
for this is that at finite k time reversal symmetry is absent 
and a non-zero value of the parameter S introduced in 
(|64|l is generally allowed. We have tuned 8 such that the 
polaronic corrections in the kinetic energy (|68a|l do not 
have the same fc-dependence as the leading term. This 
is giving rise to a a 2 / (3 term at finite k. It leads to the 
variational mass 



(m^y 1 ^ Q^E; ar (k;ao,f3 ) 



k=0 



= t 2- 



27 J 



256 V t 



(73) 



Here, the dependence of a and (3 on k 2 need not be taken 
into account because of stationarity. The coefficient of 



the second term happens to be the same as the one 
showing up in (|72|1 which is now off by about a factor of 
2 when compared to the mass shown in (|61JI . This may 
be attributed to the variational character of the state 
(|64|l . since the mass comes from a subdominant term 
proportional to k 2 . However, the correct dependence 
on J as well as the right order of magnitude are again 
obtained. We therefore conclude that the wave function 
(|FJ4"}| provides a quite accurate picture of the large po- 
laronic cloud in the limit of small J/t, especially at k = 0. 

In view of our original motivation to examine the sta- 
bility of the liquid phase of the t-J-J 1 model as x — ► 
at small J/t, it is interesting to think about the possibil- 
ity of the formation of bound hole states. It is generally 
expected that either at the critical value for the onset 
of phase separation, (J/t) c , or at an even smaller criti- 
cal value (J/t) Cl < (J/t) c bound states of two holes will 
existAi (see also Ref. 0). The existence of such bound 
states can be discussed on a qualitative level based on 
the variational spin polaron picture proposed in this sec- 
tion. To form a bound state, the single polaron wave 
functions must significantly overlap, hence the size of a 
bound state will be of order r ~ (J/t) - 2. The potential 
energy gain will be of order (J/t) s since the mean field 
term of order J in (|72|l will not be affected by pair for- 
mation. However, the kinetic energy cost of such a state 
is of order 1/r 2 ~ J/t and is dominant. We conclude 
that bound states of holes will require finite J/t of order 
1 or greater, in agreement with the picture of free sin- 
gle hole-like charge degrees of freedom established in the 
preceding sections. 



VI. CONCLUSION 

We have examined the t-J-J'-model in one dimension, 
in the regime of small x and J/t by perturbative and 
variational approaches. This parameter regime is most 
challenging to numerical methods, and earlier numerical 
studies did not allow a firm conclusion whether a phase 
separation instability and a phase of dominant singlet 
superconducting fluctuations extend down to values of 
J/t < 1 in the case a = J' / J ~ .5, where a spin gap is 
present at small doping. 

Using an approach where couplings between spin and 
charge degrees of freedom are treated as a perturbation, 
we have presented a detailed analysis of the model in 
second order perturbation theory, showing that no insta- 
bility is present at small J/t. Instead, using Luttinger 
liquid arguments and by studying the dispersion of a sin- 
gle hole immersed into the correlated spin system, we 
have demonstrated that the hole degrees of freedom pre- 
cisely behave as free spinless solitons in the dilute limit, 
despite their microscopic coupling to the non-trivial spin 
background. This behavior conforms to Luttinger liq- 
uid physics, where spin and charge are separate degrees 
of freedom, and couplings between them are regarded as 
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irrelevant in a renormalization group sense. While this 
point of view is generally accepted for one-dimensional 
systems, in microscopic one-dimensional lattice models 
it usually may be firmly demonstrated only at special 
integrable point o 2 ! 12 . The method we established in sec- 
tion IIIII provides a perturbative framework for such a 
demonstration in a non-integrable model over a range of 
parameters. Moreover, it allows the calculation of non- 
trivial quantities such as the leading corrections to the 
single hole energy and mass renormalization, which de- 
pend on non-analytic powers of J/t. The numerical cal- 
culation of the coefficients in this expansion still requires 
an exact diagonalization of a pure spin problem. We used 
these results for a comparison to a variational approach. 
Proposing a variational wave function where the hole is 
surrounded by a polaronic cloud of tightly bound pairs 
of triplet excitations we were able to confirm the pertur- 
bative results for the dependence of the single polaron 
energy and mass on J/t, as well as the order of magni- 
tude of the coefficients. In particular, the second order 
perturbative energy corrections are in close quantitative 
agreement with the variational result. Based on these 
findings, we argue that for the parameter a = .5 the on- 
set of phase separation at small doping as well as the 
formation of bound states require J/t to be at least of 
order 1. 

To conclude, the theory of the t-J- J'-model presented 
here is limited to small values of J / 1 and cannot access a 
region of dominant superconducting correlations, which 
at a = .5 might exists at moderate values of J/t and 
small doping from numerics. In a real system, additional 
interchain effects must be taken into account, that can 
either favor dimer locking, or a superconducting dimer 
liquid when dimer locking is frustrated (cf. Ref. l2ll) . 
This calls for further investigations. 
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APPENDIX A: CONVERGENCE AND 
CROSSOVER BEHAVIOR OF THE 
PERTURBATIVE EXPANSION 



Here we briefly illustrate the behavior of our expansion at 
/c'th order perturbation theory, where one will encounter 



terms analogous to (|3"Tjl : 

dhi . . . dh m / dpi... dp n 

-kf J — kf J kf "kf 

\ - J k M (hi... h m ;pi . . .p n ;oi . ..<J k -i) 
k — 1 factors 

(Al) 

The phase space consists of m-hole momenta and n- 
particle momenta. It is enough to consider the case 
m + n = k. There will be terms with fewer integrals 
also, but they are multiplied by additional powers of x 
such as pn)l. 

We focus on the regime k 2 <C J/t first. Since e(k) w 
const + tk 2 , the integrand does not significantly depend 
on the hole momenta such that each of the hole integrals 
will give rise to a factor of x. The integral over parti- 
cle momenta pt will be dominated by the region where 
all momenta are within a range of yj A/t of the Fermi 
points, where A ~ J is the spin gap. In this regime, all 
of the k — 1 factors in the denominator are dominated by 
the spin gap and are of order A. Hence we obtain the 
following estimate for the term displayed in I|A1(I 

m ~ * m (AA)" /2 ^ ~ x m (y/JJi) n J (A2) 

The leading contribution to Ek in the limit x 2 <C J /t will 
thus be a term of order 

E k ~ xU/Jft) 1 J (A3) 

Eq. (|A3(I shows that subsequent orders in perturbation 
theory are always suppressed by powers of J J/t, as we 
verified explicitly up to second order (cf. (JH3J). Note 
that relations (|A2|) and (| A3|) are valid asymptotically in 
a given limit, they do not imply the e xiste nce of a sys- 
tematic expansion in powers of x and yj J/t. Rather, the 
E k s are quite complicated functions of x and J. Rela- 
tion (|A3(I will hold until x 2 / J ~ 1/t, and upon further 
increase of this ratio a crossover will take place. We may 
however write down an asymptotic expansion in x: 

E k = L (A k (J)x + B k (J)x 2 + C k (J)x 3 + D k (J)x 4 . . .) 
E = const. + Ei + E? + . . . 

(A4) 

as we did in second order perturbation theory. Recall 
from {□), JSSJ that A 2 - (J/t) 3 / 2 , C 2 - (J/t) 1/2 while 
£?2 = 0. Formally, however, B 2 is of order J. This im- 
plies that in general yjtx 2 / J is the expansion parameter 
of the series (|A4|) . 

In the opposite limit J <C tx 2 it is easily seen from 
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(|A1|) that now Ek ~ J k holds. In this limit it is not 
necessary, though still permissible, to include the spin 
chain part H s in the unperturbed Hamiltonian Hq. In- 
stead one may apply degenerate perturbation theory in 
the spin couplings, which gives rise to an asymptotic ex- 
pansion in J: 



where H s is as defined in (J3J. Let £ q (X) denote the 
ground-state energy per site of this problem. Then it 
is easily seen from second order perturbation theory and 
the definition (|B1|I that at q = 



(B3) 



E = const + L (a(x)J + b(x)J 2 + c(x)J 3 + d(x)J 4 . . . ) 

(A5) 

This method has been applied in Ref. [l0| to calculate 
the intersection of the spin gap phase boundary with the 
x-axis. Note that in l|A5J) x need not be small, whereas 
in (| A4|> both J/t and x 2 t/J have to be small. However, 
due to the limitation J <C tx 2 , eq. (|A5|) cannot be used 
to address the nature of the phase diagram in the dilute 
hole limit. 



APPENDIX B: THE CONTINUITY OF THE 
FUNCTION F(q) 

We will now give an argument for the continuity of the 
function F(q) which leads to the crucial cancellation in 
(j36|) . This question is more subtle than it may seem, 
and the following argument would require more scrutiny 
in the gapless case a < a c . We restrict ourselves to 
the spin gapped case, as we have done throughout the 
paper. Recall the definition of F(q) from (|33J) : 



holds, where the prime denotes a derivative with respect 
to A taken at A = 0. On the other hand, at q ^ the 
same argument gives 

F{q£0) = -£>; (B4) 
Note the factor of 2 difference between (|B3|I and (|B4ll . 
Despite this apparent difference between the cases q = 
and q ^ 0, it is £'' which is discontinuous at q = 0, not 
F(q), as the following argument shows: In the vicinity of 
a site j the ground state of H q (X) will have great overlap 
with the ground state of H q= o(X(j)) as q — » 0, where 
X(j) = Acos(gj). In other words, as q — + it should be 
justified to replace the oscillating perturbation in H q (X) 
by a flat perturbation in a sufficiently large local region 
around each site j. The size of this region can still be 
chosen to be <C 1/q. One can thus argue that up to 
powers of q the ground-state energy will be given by a 
sum over local contributions £ 9= o(A(j)): 



1 1 f Ne 

£ q ^o(X) = — ^2£ q=0 (X(j)) w — / dx £ q=0 (Xcos(qx)) 



F(<l) = 



2N K 




OQ I 7-9 W) H lq K) 

E„ — E an 




1 f Ne ( 1 
— J dx [£ q=0 + X cos(qx)£' q=0 + - A 2 cos 2 (<jx)^' =0 



1 



£ q =o + -^X 2 £' q ' =Q 



(B5) 



Physically, the continuity of F(q) can be seen by in- 
terpreting F(q) as the second order energy response of 
a pure spin chain due to a periodic perturbation. More 
precisely, we consider the following auxiliary spin chain 
problem: 



H q (X) = H s + X J c os(qjhj 



(B2) 



Hs + 2 J (7 9 + 1-q) 



From JB4J, (JBSj) it then follows that 



F(q - 0) 



2 9=0 



F(0) 



(B6) 



Note that the local point of view taken here is better 
justified in the gapped case, where any local perturbation 
decays exponentially in space. 



1 J. Voit, Rep. Prog. Phys. 58, 977 (1995). 

2 F. D. M. Haldane, Phys. Rev. Lett. 45, 1358 (1980). 

3 M. Imada, Prog. Theor. Phys. Supp. 113, 203 (1993). 

4 M. Nakamura, J. Phys. Soc. Jpn. 67, 717 (1998). 

5 M. Ogata, M. U. Luchini, S. Sorella, and F. F. Assaad, 
Phys. Rev. Lett. 66, 2388 (1991). 

6 E. Dagotto and T. M. Rice, Science 271, 618 (1996). 

7 M. Imada, Phys. Rev. B 48, 550 (1993). 

8 A. Tanaka and X. Hu, Phys. Rev. Lett. 88, 127004 (2002). 

9 P. W. Anderson, Science 235, 1196 (1987). 



M. Ogata, M. U. Luchini, and T. M. Rice, Phys. Rev. B 
44, 12083 (1991). 

1 F. D. M. Haldane, private communication. 

2 H. J. Schulz, Phys. Rev. Lett. 64, 2831 (1990). 

3 J. Solyom, Adv. Phys. 28, 201 (1979). 

4 R. Jullien and F. D. M. Haldane, Bull. Am. Phys. Soc. 28, 
344 (1983). 

5 F. D. M. Haldane, Phys. Rev. B 25, 4925 (1982). 

6 T. Xiang and N. dAmbrumenil, Phys. Rev. B 45, 8150 
(1992). 



14 



W. Kohn, Phys. Rev. 133, A171 (1964). 
A. Seidel, Ph.D. thesis, Massachusetts Institute of Tech- 
nology (2003). 

C. K. Majumdar and D. K. Gosh, J. Math. Phys. 10, 1388 
(1969). 

W. J. Caspers and W. Magnus, Phys. Lett. 88A, 103 
(1982). 

A. Seidel, C. A. Marianetti, F. C. Chou, G. Ceder, and 
P. A. Lee, Phys. Rev. B 67, 020405(R) (2003). 



We thank F.D.M. Haldane for pointing out to us this pos- 
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Note that for small doping and J/t the (PS) region consists 
of an insulating x = phase (no holes) and a liquid phase 
with x = x c , where x c is the critical doping at which phase 
separation just occurs. Since x c is small in the region of 
interest, it is justified to consider the holes as dilute even 
in the (PS) regime. 



